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ABSTRACT 



We present infrared aperture synthesis maps produced with the upgraded 
IOTA interferometer. Michelson interferograms on the close binary system 
Capella (a Aur) were obtained in the if -band between 2002 November 12 and 16 
using the IONIC3 beam combiner. With baselines of 15 m < B < 38 m, we were 
able to determine the relative position of the binary components with milliarc- 
second (mas) precision and to track their movement along the ~ 14° arc covered 
by our observation run. We briefly describe the algorithms used for visibility 
and closure phase estimation. Three different Hybrid Mapping and Bispectrum 
Fitting techniques were implemented within one software framework and used 
to reconstruct the source brightness distribution. By dividing our data into sub- 
sets, the system could be mapped at three epochs, revealing the motion of the 
stars. The precise position of the binary components was also determined with 
model fits, which in addition revealed lAal^Ab = 1.49 ±0.10 and apparent stellar 
uniform-disk (UD) diameters of Qah = 8.9 ± 0.6 mas and ©a& = 5.8 ± 0.8 mas. 

To improve the u, f-plane coverage, we compensated this orbital motion by ap- 
plying a rotation-compensating coordinate transformation. The resulting model- 
independent map with a beam size of (5.4 x 2.6) mas allows the resolution of the 
stellar surfaces of the Capella giants themselves. 

Subject headings: techniques: interferometric - methods: data analysis (wavelet) - 
binaries: close - spectroscopic - stars: binaries - close - spectroscopic - stars: individual 
(Capella) 
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1. Introduction 

Aperture synthesis images, obtained using an interferometer with multiple baselines 
and closure phases, provide a much more assumption-free way to interpret interferometric 
data astrophysically than the simple fitting of models to single-baseline fringe amplitudes 
does. Therefore, imaging capabilities will be of special importance for the upcoming 
generation of ground-based and space interferometers. In this paper, we present the results 
of an initial attempt to use the Infrared Optical Telescope Array (IOTA) for imaging in the 
near infrared. Since this was to be the first effort to construct an image at IOTA, we decided 
that it was important to observe an object with a well known brightness structure in order 
to demonstrate the instrument's capabilities. Therefore, we chose to observe the strongly 
resolved spectroscopic binary Capella (a Aurigae) so that our maps could be compared to 
the results of earlier studies. 

The first spectroscopic measurements by Campbell (1899) revealed the components 
Capella Aa and Ab (classified today as G8 III and Gl III (Johnson et al. 2002)). The system 
also became the first binary to be separated by an optical interferometer (Anderson 1920). 
The orbit was determined with high precision in 1994 using the Mark III interferometer 
by Hummel et al. (1994). In 1996 Baldwin et al. presented the first optical (A = 830 
nm) aperture synthesis map of Capella using the Cambridge Optical Aperture Synthesis 
Telescope (COAST) with a longest baseline of B = 6.1m. Two years later Young used the 
same instrument to map Capella in the near infrared (A = 1.3 /iin; Young 1999). 

We observed Capella on five nights in November 2002 in the if-band (A = 1.65/iin; 
A A = 0.3/iin). In this paper, we review the data reduction procedure used to determine 
visibilities and closure phases from the raw interferometer data, and we present model 
fits and images derived from the visibility and closure phase data. By investigating the 



-4- 



well-studied Capella system, we have been able to calibrate the closure phase sign and to 
determine the overall behaviour of this recently upgraded 3-telescope interferometer. One 
of the mapping procedures described within this paper has also been used to generate the 
maps presented in Monnier et al. (2004) depicting the binary A Virginis. 

2. Observations 

The IOTA was constructed jointly by the Smithsonian Astrophysical Observatory 
(SAO), Harvard University, the University of Massachusetts (UMass), the University of 
Wyoming, and MIT/Lincoln Laboratory (Traub et al. 2003). This 3-telescope interferometer 
located at Fred Whipple Observatory atop Mount Hopkins, Arizona, demonstrated first 
closure phases on 2002 February 25. IOTA's telescopes can be moved on an L-shaped track 
and are mounted on stations si5m and « 7m apart. With track lengths of 15 m (south-east 
direction) and 35m (north-east direction), baselines up to 38m can be formed. The collecting 
optics consist of f/2.5 45 cm Cassegrain primary mirrors, which are fed by siderostats. 
The atmospherically induced motion of the image is compensated by tip-tilt servo systems 
mounted behind the telescopes. Passing various mirrors and path-compensating delay lines, 
the beams are deflected into the laboratory, where the infrared component of the beam is 
coupled onto the IONIC3 pair-wise beam combiner (Berger et al. 2003). This combines the 
beams symmetrically with a ratio of 50:50. 

For each baseline the beam combination produces two complementary outputs which 
are simultaneously acquired by a PICNIC-camera (Pedretti et al. 2004). Using piezo 
scanners installed in two optical path delays, the optical path delay (OPD) can be 
modulated to scan temporally through the interference fringe pattern. 

Each observation on the target object is completed by the acquisition of four calibration 
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Table 1. Observation Log 



Date [UT] 


Start time [UT] 


No. exposures^ 


Calibrators 2 


11/12/02 


07:43 


20 


a Cas, Aur 




10:11 


27 


Aur 




12:00 


16 


Aur 


11/13/02 


07:19 


26 


8 Aur, Aur 




09:09 


20 


Aur, 8 Aur 




10:45 


24 


8 Aur, Aur 


11/14/02 


05:51 


25 


k Per, Aur 




08:24 


21 


5 Aur, Aur 




10:25 


39 


8 Aur, Aur 




12:06 


19 


Aur, 8 Aur 


11/15/02 


07:08 


25 


a Cas, <5 Aur 




09:28 


53 


<5 Aur 




11:34 


18 


8 Aur 




12:50 


6 


5 Aur 


11/16/02 


09:12 


45 


8 Aur, /3 Aur 




12:22 


17 


8 Aur 



Interferometer Configuration used: 
2002 Nov 12-14: A35, B15, CO 
2002 Nov 15-16: A15, B15, CO 
where A, B, C denotes the individual telescopes and the number of their position 
on the north-east (for A and C) or south-east baseline (for B) in meters. 

x Each exposure consists of 500 individual scans. 

2 Compared to other calibrator stars, a Cas is clearly resolved. In 
order to still use it as a calibrator we applied a UD correction based 
on the reference diameter © a c as = 6.25 mas given by Richichi and 
Pcrchcron (2002). Also the well-known binary Aur with a semi- 
major axis of 3.3 mas appears marginally resolved to our interfer- 
ometer. We used the orbital elements and UD diameters derived 
from measurements with the Mark III interferometer (Hummel et al. 
1995) to apply a correction to our measured visibilities and closure 
phases before calculating the transfer function. The short orbital pe- 
riod (3.96 d) requires to compute this correction for each of our mea- 
surements on Aur separately. 
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files. One of these files measures the camera background signal with the light from all 
telescopes shuttered out. For the other three files, the light of two telescopes is shuttered 
out alternately. This allows us to measure the instrument transmission by determining the 
coefficients of the transfer matrix (k matrix, Coude Du Foresto et al. 1997). 

The observations presented in this paper were obtained at the IOTA between UT 
2002 November 12 and 16. During this time we made observations in two three-telescope 
configurations in order to obtain a reasonable sampling of the u, v plane for the Capella 
image. A log of the observations is presented in Table 1. 

3. Data Reduction 

(see Figure attached) 

Fig. 1. — Representative corrected interferograms, all taken 2002Novl3 on the BC baseline. 
The selected scans represent typical conditions with low noise (a Cas), high noise (5 Aur), 
and interfering resonance (a Lyn). The wavelet spectral density power spectra shown in the 
lower panels are calculated from the corresponding raw interferograms. On the ordinate of 
the CWT plots, the scale quantity is given which is equivalent to the period and inversely 
proportional to the frequency of the corresponding rescaled wavelet. The contours demon- 
strate how our departitioning algorithm removes regions not connected to the area with the 
highest peak and separates the fringe peak from the piston and resonances for the visibility 
estimation. 

IOTA produces raw fringes by scanning the optical path delay through the position 
of the white light fringe in the instrument. The fringes which are observed must then 
be converted into visibilities and closure phases for imaging and model fitting. We have 
reduced the raw fringe data in a sequence of steps. 



-7- 



An important first step is to locate the fringes within the raw data scans in order to 
determine the optical path delay in the IOTA system. If fringes are found, the next step 
is to calibrate photometric fluctuations due to atmospheric seeing. By using the transfer 
matrix, the coupling of the stars to the optical inputs of the beam combiner and the internal 
coupling of the beams within the combiner (Coude Du Foresto et al. 1997) are taken into 
accont as well. 

In order to form a photometrically calibrated scan, we make use of the fact that for 
the 3-beam pairwise combination scheme used the photometric information can be obtained 
for each measurement using linear combinations of the interferometric signals (similar to 
the procedure used by Monnier 2001). This requires a matrix inversion, which might be 
sensitive to numerical instabilities, especially in the case of low signal-to-noise (S/N). For 
this reason, we smooth the scan by convolving it with a Gaussian before solving for the 
photometry for each individual pixel. The photometric calibration is performed for the two 
complementary signals of each baseline separately. Remaining residuals are removed by 
subtracting the two complementary outputs for each baseline. 

The upper panels of Figure 1 illustrate how the data might look for different conditions 
after these basic data reduction steps. 

Each observation at IOTA consists of measurements of hundreds of fringes obtained 
over a time period of a few minutes. Typically, we have estimated visibilities and closure 
phases for each individual realization of a fringe and then averaged the results for many 
fringes to obtain the final result. 
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3.1. Visibility Estimation 

The scans through the white light fringe follow a particular functional form. Assuming 
a bandwidth filter with idealized rectangular transmission, the instrumental response 
expected for a point source is given by the product of a sinc-envelope, with width inversely 
proportional to the bandwidth of the input signal, and a cosine function representing the 
interference fringes. The amplitude of the fringe envelope reflects the raw visibility. However, 
this measurement must be calibrated with observations of unresolved reference stars or 
stars of known diameter in order to remove instrumental and atmospheric effects. Typically, 
each observation of a target source is bracketed by observations of unresolved stars to 
provide a calibration of the instrument. In this work, all target observations are enclosed 
by at least one such calibration measurement on each side of the target observation. The 
results for the calibration stars are interpolated linearly to calibrate the raw visibilities for 
the interferometric efficiency of the system, yielding the calibrated target visibility. 

During our observation run we noticed a minor technical problem which became a 
dominating issue for data reduction: the scanning piezo mirrors showed a mechanical 
resonance (see 3rd panel of Figure 1), resulting in multiple additional peaks within the 
power spectrum of the acquired scans. The data reduction procedure most commonly used 
by the interferometry community (e.g. Baldwin et al. 1996) estimates the visibility by 
integrating over the fringe peak within interferogram power spectra after averaging and 
background subtraction. To improve the visibility estimation for our resonance-influenced 
data, we implemented two alternative methods to estimate the visibility amplitude. 

Since the interfering resonance only affects the internal structure of the fringes, we 
implemented a Fringe Envelope Fitting (FEF) algorithm which filters the peaks of the 
fringe packet. The position of these local minima and maxima is then fitted to the analytic 
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sinc-envelope function with the envelope amplitude, the width, and the position as fitting 
parameters. For low S/N this algorithm tends to overestimate the fringe amplitude. 
Analyzing the resulting bias by simulating fringes with additive Gaussian noise, we found 
a bias on the fitted fringe amplitude which depends exponentially on the noise. Based on 
these simulations we were able to correct the fitted visibilities down to S/N I. 

Motivated by Segransan et al. (2003), we also implemented a data reduction procedure 
based on the Continous Wavelet Transform (CWT). In multiscale analysis, the time signal 
is decomposed into the two-dimensional time-frequency space by measuring the response of 
a translated and scaled mot her- wavelet function to the measured signal. Compared to the 
Fourier transform (FT), which is based on periodical sinusoidal waves, the use of a localized 
mother wavelet allows the precise localization of the power within the signal (Torrence and 
Compo 1998). For the mother wavelet, we used the Morlet function, which is given by a 
sinusoidal multiplied with a Gaussian and therefore similar to the analytic fringe function 
as noted above. The absolute square of the complex wavelet function defines the wavelet 
spectral density spectrum. 

From this density spectrum, we remove all features below a particular significance level 
(e.g. 35% scaled to the peak). Resonance (see Fig. 1) or an uncompensated change in piston 
during the fringe scan might result in several non-connected areas. Therefore, we run a 
filling algorithm which removes all regions which are not connected to the area containing 
the main signal. The integral over the remaining region is used as an estimate of the fringe 
power. 

The measured width of the remaining region in the OPD and scale domain is then 
compared to the expected width and used as a selection criteria to reject scans with a high 
piston. Even for good scans, the CWT is contaminated by additive noise. To determine the 
level of this noise, we integrate in the CWT along the OPD axis towards both directions, 
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starting 5/xm from the border of the fringe region. The extend of this noise region along 
the scale axis is given by the extend of the fringe peak region. This averaged background 
noise level is subtracted from the fringe power. Noticing that the CWT is a non-linear 
transformation, we confirmed the proper behavior of our algorithm by performing extensive 
simulations. 

The visibilities obtained with the FEF and CWT algorithm are in good agreement, 
even if the CWT algorithm seems to be superior for cases with low S/N. The accuracy 
of the measured visibilities is limited primarily by calibration. For error estimation we 
include, beside the statistical errors for each data set, a calibration error for each night 
which was estimated by the scattering of the calibrator visibilities during that particular 
night. For good seeing conditions this calibration error AV^ aHb is typically around ±1% 
(2002 Nov 12-15), while for worse conditions it may rise to ±2% (2002 Nov 16). 

3.2. Closure Phase Estimation 

Since we are observing on the three baselines AB, BC, AC simultaneously, the fringe 
phases 0ab, <Pbc, <Pac can be measured. Simultaneous measurement of three phases in turn 
allows a new quantity, the closure phase $ (CP, Jennison 1958), to be estimated. The CP 
has enormous value since it is unaffected by atmospheric propagation effects, and through 
measurements of CP and visibility amplitudes, it is possible to create aperture synthesis 
images. In the following we briefly describe our methods for estimating this essential 
interferometric observable. 

The estimation of CP begins with the location of the fringes for each baseline within the 
raw scans. The phase of each fringe must be determined within the same temporal window 
within the scan in order to assure that all phases are measured at the same time and within 
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the coherence time of the atmosphere. Thus, we locate the median delay for the centers of 
the three fringe packets and determine the phase of each baseline within a short window of 
width 3/iin centered on that OPD. The FT of the windowed fringe results in a strong peak 
at the frequency of the fringe. In the IOTA three-baseline system, the frequencies of these 
peaks for the three fringes obey a "closure" relationship such that uab + ^bc — v ac = 0, and 
if this relationship is not obeyed by the data, then there is a good chance that one or more 
of the peaks is spurious. For these reasons, we reject scans which do not obey this frequency 
closure relation. For scans that do pass this test, we may then use the phases of the FT's 
for the three fringes to compute the CP: $ = 4>ab + 4>bc — 4>ac- Unresolved calibrators, for 
which the CP is known to be zero, are used to define the instrumental CP-offset. 

Using the above algorithm, we have estimated closure phases for all of the observations 
obtained during our run. From the scattering of the calibrator CP over the individual nights 
we include calibration errors A<3> ca / ife between ss 0.2° and 1.1°. Under good seeing conditions, 
the IOTA/IONIC3 instrument is able to measure CP with high precision with systematic 
errors of the order of 1.0° or less. As a conservative value we assume A$ SJ/S = 2.0°. 

3.3. Investigating the Effect of Bandwidth Smearing 

After presenting the details of our procedures for visibility and CP estimation, we want 
to investigate a potential influence of the bandwidth smearing effect on our data. 

For extended sources, this effect will reduce the complex degree of coherence at those 
points of the brightness distribution that are separated in delay from zero optical path 
delay r by amounts comparable to the width of the fringe packet. As a consequence, the 
detailed instrumental response will be more complicated than the point-source response. 
Thus, estimation of visibilities and closure phases using the point source response as a 
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template could lead to some systematic errors. 

For a rectangular spectral bandpass profile, the complex degree of coherence c is given 



by 



_ sin(7rrAA/A 2 ) 
{ > tttAA/A 2 • 1 } 



To estimate the loss of coherence for a particular projected baseline B±, the coherent field 
of view Aa = r/B± can be defined as the angular separation for which the fringe of an 
off-axis point source will be centered right on the first zero point of the fringe envelope of 
the on-axis point source. In this case, just half of the fringe power from both sources will 
overlap. Using Equation 1, we obtain Aa = A 2 /(_Bj_AA) 49 mas for our longest baseline 
used (AB on the nights of UT2002Novl2-14). The comparison with the separation of the 
Capella stars (s ~ 46 mas) shows that we just reach the limit where bandwidth smearing 
might play a role. 

Since in the regime of marginal bandwidth smearing the contrast decreases with 
distance from zero OPD in a well-defined way, a correction can be applied. In the given case 
we correct the visibility Vud(Q) of a uniformly bright disk of diameter 9. Assuming totally 
coherent combination, Vjjd can be written as 

Vud{0) ~ kQB ± ' (2) 
where J\ denotes the Bessel function of first kind and first order. 

To determine the delay- affected visibility , we define 

V^T(r) = c{r)V UD (3) 

In the case of a binary star with two UDs of diameter Qao, and &Ab the resulting 
interferogram is given by the composition of two fringes with different delays r Aa and r Ab . 
Since the delay at which the combined fringe is evaluated depends on the brightness ratio 
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of the individual fringes, it is appropriate to evaluate the two correction factors relative to 
an effective delay r e g, given by the center of intensity of the combined fringe packet: 

I AaTAa + I Ab^Ab /-x 

Teff = I Aa+ I Ab (4) 

v UD,Aa = c(r eS -r Aa )V UD (e Aa ) (5) 
VJZm = c(r eS -r Ab )V UD (e Ab ). (6) 

The fact, that our procedure for closure phase estimation does not measure the phases 
over the whole fringe packet could be of advantage in the regime of marginal bandwidth 
smearing. As described in chapter 3.2, we measure the phases within a narrow temporal 
window around the center of the composite fringe where the interferograms overlap properly. 



4. Model Fitting 

Making use of the expression for the visibility of a delay-affected UD, the complex 
visibility V of a close binary system with a position vector x can be written as 

v /^\ _ I Aa v UD,Aa + 1 AbVup )A b e /„s 

V \ X ) — J , J 5 {') 

i A a + iAb 

— * 

with B = (u,v) for the baseline vector and the wave number k = 2n/X. In strictly formal 
notation, an additional overall phase factor appears on the right-hand side of Equation 7. 
Since this phase factor merely moves the fringe pattern underneath the envelope it neither 
effects the visibility nor the CP and can be neglected. 

We performed least square model fits, equally weighting N v individual visibility and 
CP measurements with 



X 2 = X 2 V 2+xl (8) 

Nv / T/2 _ T/2 \ 2 



i=l 
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xl = E( 

i=l x 



2 



(10) 



with the model visibilities V mo d e i and closure phases Q mo dei calculated from the complex 
visibility: 



The errors for the individual measurements are ayi and <r$. The model for Capella included 
the relative positions of the stars, the stellar diameters, and the ratio of the total fluxes 
from the two stars. 

Given our limited dataset, the Capella model contains parameters with significant 
correlations, which makes it difficult to extract all parameters simultaneously. Therefore, 
we proceeded to estimate the model parameters in an iterative manner, starting without 
delay-compensated visibilities (using V UD instead of Vjf^ in Equation 7). In the first step, 
the stars were assumed to be point sources and the fit solved for the relative positions and 
fluxes. The next step was to fit the residuals to this model fit and determine the stellar 
diameters. The process continued iteratively, using the derived diameters to find relative 
fluxes and positions followed by a new estimation of the diameters. This process converged 
typically within a few iterations. Finally, we applied the delay-compensation to obtain the 
fit results for each individual night as given in Table 2 and Fig. 2. 

The fit results for each individual night are given in Table 2 and Figures 2 and 3. In 
some cases, the reduced x 2 values exceed unity, which might be a result of bad seeing, a 
potentially uncorrected bias due to the effect of bandwidth smearing, the piezo scanner 
resonances, or other undiscovered systematic errors. However, the measurements are clearly 
in good agreement with the high precision reference orbit by Hummel et al. (1994) obtained 
with two telescopes on the Mark III interferometer. 
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Table 2. Fitting Results 



Data from Night [UT] 


Data Points 


X 2 /DOF 


Fitted 1 


Fitted Diameter 1 


Fitted Position 2 


Ref. Pos. 2 ' 3 


Date MJD=JD-2452500 


N v N 9 


Xv x% 
N v Nq, 


lAb 


B Aa [mas] B A b [mas] 


dRA [mas] dDEC [mas] 


dRA [mas] dDEC [mas] 



11/12/02 


90.821 . 


. 91.046 


165 


54 


0.76 


1.15 


i 94+ - 36 


8.6: 


-1.0 

-0.8 




-9.28 ± 1.94 


44.79 ± 1.25 


-8.93 


44.70 


11/13/02 


91.805 . 


. 91.978 


173 


59 


1.11 


2.49 


i r-,+0.29 
J - ol -0.31 


9.6: 


-1.4 
-1.2 


5 2+ 3 - 4 


-12.14 ± 1.68 


43.33 ±0.89 


-11.67 


43.25 


11/14/02 


92.744 . 


. 93.026 


315 


104 


1.05 


3.75 


1 42+ - 25 
^ -0.20 


8.7: 


-1.1 

-1.0 


6-o±l:S 


-14.29 ± 1.24 


41.32 ±0.71 


-14.47 


41.60 


11/15/02 


93.797 . 


. 94.042 


363 


105 


0.29 


0.75 


i crn + 0.20 
i - ou -0.10 








-17.90 ±3.90 


39.27 ±4.43 


-17.33 


39.72 


11/16/02 


94.883 . 


. 95.034 


183 


58 


0.56 


0.24 


-, co+0.29 
1 - b8 -0.26 








-20.24 ± 1.40 


37.68 ± 2.40 


-20.13 


37.68 



1 To compensate the perceivable intra-night motion, the same rotation-compensating coordinate-transformation as described in chapter 5.3 was applied separately to 
the data of each night before fitting the intensity ratio and the diameter. 

2 Relative positions are measured from the infrared brighter to the infrared fainter component. 

3 Rcference orbit by Hummel ct al. (1994). 
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Fig. 2. — Best-fit models for Capella assuming uniform disks with parameters as listed in 
Table 2. For all plots, the error bars shown include statistical, calibration, and systematic 
errors as discussed in chapters 3.1 and 3.2 and the CP is given in radians. The abscissa 
depicts MJD=JD-2452500. 
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Fig. 3. — Best-fit models for Capella assuming uniform disks with parameters as listed in 
Table 2. For all plots, the error bars shown include statistical, calibration, and systematic 
errors as discussed in chapters 3.1 and 3.2 and the CP is given in radians. The abscissa 
depicts MJD=JD-2452500. 
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Weighted averaging give for the diameters of the Capella giants, &Aa = 8.9 ± 0.6 mas 
and 0^6 = 5.8 ± 0.8 mas and the intensity ratio iAa/^Ab = 1-49 ± 0.10. Using the orbital 
parallax of 75.0±0.57mas (Pourbaix 2000) we calculate physical radii of R Aa = 12.6±0.9 R f 
and R Ab = 8.3 ± 1.1 R Q . 

Stellar radii measurements are of special importance as the most direct way to 
determine stellar effective temperatures. The stellar radii, derived from our H-ba.nd data, 
are in good agreement with the measurements at visible wavelengths by Hummel et al. 
(1994). As a result of our longer effective wavelength and the smaller amount of data, we 
cannot exceed the precision of their extensive study. Therefore, we skip re-calculating the 
effective temperatures with our values. 



5. Imaging 

5.1. Imaging Algorithms used 

In the exploratory spirit of the present Capella imaging work, we have investigated 
three different imaging algorithms, as explained here. 

The radio astronomy community has developed various algorithms to reconstruct phase 
information from CP, most notably conventional hybrid mapping (CHM, e.g. Cornwell 
and Wilkinson 1981) and difference mapping (DFM) by Pearson et al. (1994). In both 
algorithms, the individual phases are recovered by self-calibration (Pearson and Readhead 
1984). The necessary deconvolutions can be performed with the classical CLEAN algorithm 
(e.g. Hogbom 1974). 

We have also tested the "building-block mapping" (BBM) algorithm presented 
by Hofmann and Weigelt (1993), which was originally developed to reconstruct images 
from speckle-bispectrum data. Here, the image is build up by adding iteratively point-like 




Fig. 4. — Image reconstructions from subsets of our complete data set revealing the or- 
bital motion of the stars between three epochs. To compensate the marginal motion within 
the subsets, the u, f-plane was rotated synchronously to the reference position of the 
components at the intermediary modified julian dates MJD=91.437, 93.907 and 94.262 
(MJD=JD— 2452500). The crosses within the maps indicate the "centers of light" of the 
stars within the other images. Whereas the positions fit the expectations accurately (< 1 
mas), the intensity ratio between the two stars was not always obtained properly due to the 
poor sampling (the coverage is shown below each map). The contours show 10% intervals 

snnlprl tn tViP npnlr 
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Hybrid Mapping Difference Mapping 




dRA [mas] u [10 6 A] 

Fig. 5. — Aperture synthesis maps of Capella generated with data from 2002 Nov 12-16. Maps 
generated with the CHM (upper left), DFM (upper right), and BBM (lower left) algorithm 
are shown together with the u,v -plane coverage (lower right). Each map shows the surface 
brightness plotted with 5% interval contours (scaled to the peak intensity). To compensate 
for the motion of the stars over the observed interval (~ 14°), the u,v -plane was moved 
synchronously to the orbital motion. The maps are oriented north up, east left. The restoring 
beam (5.4 x 2.6 mas) is shown in the upper-left corner of each map as inset. 
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components (building blocks) to an image in order to minimize the least square distance 
between the measured bispectrum and the bispectrum of the reconstructed image. Since 
the computation of this distance function for each potential position of the next building 
block (the whole image space) is computationally too expensive, a linear approximation of 
the least square distance is applied. The main assumption of this approximation is that the 
added component changes the bispectrum of the reconstructed image only insignificantly. 
For details of this method, see Hofmann and Weigelt (1993). 

In order to allow a comparison of the results obtained with the various methods, we 
implemented all algorithms within one software framework. This allows us to start all image 
reconstruction methods with precisely the same initial conditions. 

5.2. Mapping the orbital motion 

To demonstrate that the data and mapping techniques were sensitive to the movement 
of the stars during the observation period, we subdivided our data into three subsets 
(2002 Nov 12-13, 2002 Nov 14, 2002 Nov 15-16), each one still containing enough data for 
mapping. Since we observed during the first two epochs with the telescope configuration 
including the longest baseline (and hereby missing lower frequencies), the convergence of 
the BBM maps had to be supported by limiting the reconstruction area. Within the maps 
generated (Fig. 4) it is perceivable that the orbital motion can clearly be traced accurately. 
Also the intensity of the stars (from left to right: iAa/^Ab ~ 1-6, 1-7, 1.6) can be estimated 
reliably. 
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5.3. Compensating the orbital motion 

The key goal of our project was to produce aperture synthesis images of Capella in order 
to demonstrate the capabilities of IOTA/IONIC3. However, observations were obtained 
over five nights, during which time the relative positions of the binary pair changed by 
many milliarcseconds. In order to make a single map with data from all nights, and thereby 
improve the u, v coverage available for the map, we applied a coodinate transformation 
to compensate for the orbital motion of the binary components over the observed time 
interval. Thus, the fringe spacing is rescaled to account for the changing separation of the 
stars, and the orientation of the fringes on the sky is reoriented to account for the changing 
position angle of the binary pair. Each u, v point was transformed in this manner to a new 
coordinate system corresponding to the geometry of the stars at the time of the middle of 
our observation run. This way, all observations may be combined to produce a single map 
with improved u, v coverage. 

Of course, the above procedure is only strictly justified for special circumstances. In a 
binary system, the stars would have to rotate synchronously and any brightness structure 
on each star would have to be stationary during the time of observation. In the Capella 
system, we note that this procedure might be justified for the cooler Aa component, which 
has been found to rotate synchronously to the orbital motion (Strassmeier et al. 2001). The 
infrared-fainter Ab component rotates asynchronously, hence it would not be possible to 
map structure onto this more active star. 

The resulting CHM, DFM and BBM maps are presented in Fig. 5. For the given, 
sufficient simple brightness source distribution, all algorithms converged without any 
concrete trial model assumptions or other specific pre-suppositions. As an initial phase 
estimate we set two of the three phases equal to zero, thus the third phase was given by 
the closure phase relation. The intensity ratios iAa/^Ab measured within the maps are 1.9 
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(CHM), 1.6 (DFM), and 1.8 (BBM). 

A quantitative comparison between the BBM map in Fig. 5 and the maps constructed 
from the data subsets as described in chapter 5.2 reveals that the rotation compensating 
coordinate transformation clearly reduced the noise content within the maps and also 
improved the convergence behaviour. 

Since the stellar surfaces within the maps are partially resolved, one can measure 
the stellar diameters. This provides values of 0^ a ^ = 8.3 ± 1.6 and @^ a ^ = 6.8 ± 1.2 
(measured within the unconvolved BBM map), which is in good agreement with our model 
fitting results. 

6. Conclusions 

By observing the well-studied Capella system we were able to demonstrate the overall 
performance and imaging capabilities of the IOTA/IONIC3 instrument. IOTA was able to 
successfully measure and calibrate visibility amplitudes and closure phases, even with a 
temporary but significant instrumental problem that affected the raw data adversely during 
this run. The visibility amplitudes and closure phases were then used with three algorithms 
to reconstruct reliable images of this well known binary system. Thus, the primary goal of 
this experiment has been achieved, and the IOTA facility has been shown to be a powerful 
tool for high-resolution infrared imaging. 
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